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The tilted axis cranking model is used in combination with the random phase approximation 
and particle number projection to analyze the influence of dynamical pair correlations in the high- if 
bands of 178 W and their effect on relative energy and angular momentum. The calculations show 
the importance of dynamical pair correlations to describe the experiment as well as advantages and 
problems with the different models in the superfluid and normal state regions. 
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The transition of a nucleus from the superfluid to the normal state at high angular momentum is an interesting 
problem that is studied by means of modern 7-detector arrays. Contrary to the analogous transition in solids, in the 
finite nuclear system there is no sharp phase change but an extended transition region with in which pairing effects 
disappear. The most rapid attenuation of pair correlations is caused by quasiparticle excitations i.e. breaking of 
pairs. This case is realized in deformed nuclei when a large fraction of the angular momentum is generated along the 
symmetry axis of the nucleus. These states may appear as high- if isomers near the yrast line. Hence, the experimental 
data on high- if isomers and the rotational structures built on those states contain valuable information about the 
pair correlations and how they are influenced by the rotation. 

The theoretical analysis of high- if band structures presented in this paper is based on the tilted axis cranking 
(TAC) model, which is a mean field approach for describing both the rotation and pair correlations in the framework 
f*"** , of the Hartree-Fock-Bogoliubov (HFB) theory, fn fact, the TAC model has already been applied to the high- if 
multi-quasiparticle bands in 178 ' 179 W There, only the effects of static pair field have been considered. The central 
aim of our present investigation is studying the role of dynamical pair correlations. For this purpose we apply the 
random phase approximation (RPA) to the pairing interaction, which includes the fluctuations of the pair field. The 
fluctuations are particularly important when the static pair field has collapsed. Their relevance is suggested by the 
results of |2j|j| for 178 < 179 W as well as in the earlier investigations of the pair correlations of other high- if band heads 
states 0. 

The combination of the HFB theory with RPA does not provide an reliable description in the region where the 
pairing gap disappears f|. There the particle number projection method (PNP) works better PH]. Therefore, both 
^ , methods, RPA and PNP, are considered and compared with each other. 

The paper is organized as follows. In section |l| we develop TAC versions which include pairing- RPA or PNP. Details 
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of the calculations are given in section III. The results of the calculations of high- if bands in 78 W are presented in 



section IV and compared with the experiment 



II. THE MODEL 



A. Mean field 



Our investigations are based on the single particle Routhian |jj of the TAC model complemented by a monopole 
pair interaction term, 

H' = t + V{e v ) - w(sinz9j! + costfj 3 ) - GP^P - \N, (I) 

where P^ denotes the monopole pair field operator Jsj and N is the particle number operator. For simplicity the terms 
are written only for one kind of particles. 

Replacing in 7i' the pairing two-body term by the pair potential one obtains the quasiparticle Routhian 

ti = t + V (e u ) - w(sin §j x + cos <&j 3 ) - A(P* + P) - \N. (2) 

The diagonalization of h! (for details cf. ||) provides the quasiparticle energies e'^ui) as well as the HFB amplitudes 
(it, v), required later on for the RPA and PNP calculations. In order to define a quasiparticle state | > for a specific 
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configuration the occupation numbers need to be chosen. The self-consistent treatment of h' implies determining the 
pair gap A from the HFB gap equation 

G < P f >= A (3) 
and fixing the chemical potential A by the particle number condition 

< N >= N. (4) 



These conditions remain valid for the RPA but they will be modified for the PNP approach (cf. section II C). The 
HFB part of the total Routhian becomes 

Rhfb =< W > +XN. (5) 

In order to determine the deformation parameters e v we have constructed the Total Routhian Surface (TRS) using 
the standard Strutinsky renormalization procedure as described e.g. in |||]]. As snown m JltJ the essential shell- 
correction part Egtruti^u) of the TRS can be calculated for a non-rotating and unpaired ground state. The contribution 
containing the dependence on rotation and pair field enters via the quasiparticle Routhian, eq. (||). Hence, the total 
Routhian can be written as 

R(u, e v , A, A, d) = E strut (e v ) + R H fb(u, e„, A, A, ■&). (6) 



B. Pairing- RPA 

So far only the mean field part of the pairing energy has been taken into account. Now we include the pair 
correlation energy, which comes from quantum vibrations around the mean field minimum. The RPA treatment of 
the pair interaction gives the following expression for the correlation energy || 



E RPA = - 
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+ E ex (7) 



where £l v are the RPA frequencies for the pair vibrations, = + ej (/i = i < j) are the two-quasiparticle energies 
and E ex is the so-called boson exchange term. The total Routhian is the obtained from eqs. (^|,0) as 

R rpa {lo, £„, A, A, 0) = R(u, e„, A, A, d) + E* PA (u, e v , A, A, 0). (8) 

Due to the large number of RPA roots in deformed nuclei, it is practically impossible to evaluate the sum in eq. (0) 
directly HI). However, it can be calculated by means of the integration method developed recently (ll]], which is 
especially simple for pairing interaction. It uses the RPA response function F(Ci) (cf. HH), which reads explicitly for 



our case: 



where 



S ^ = ^§ Yl i ( Uk * U kj ± VkjVkv) ~ (ukjUki ± V kl V- kj ) } (10) 

and u, v are the HFB amplitudes. The zeros of the function F(f2„) determine the the RPA frequencies, Q, v . Continuing 
the variable Q into the complex plane z, one defines the spectral function ■prjr- It has first order poles at £l u and E^ 
and it is analytical for all other complex values of z. According to Cauchy's theorem the following integral relation is 
obtained: 

1 i^^n*) 
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where g{z) is a arbitrary complex function which is analytical within the region enclosed by the integration path C. 
The roots i7„ and the poles of F(z) lie in the same region. By choosing g(z) = z the RPA correlation energy, 
eq. (0) becomes 

E ^=hh*w> +E -' (12) 

where the integration path C goes around the right half of the complex plane. The exchange term is given by 

E <* = % E ( 13 ) 
n 

as in || . The integral in eq. (|ll]) is independent of the path C as long as all poles in the positive plane are enclosed 
(see Fig. [[]). Therefore C can always be chosen in such a way that the spectral function becomes smooth and the 
integration numerically stable, such that a small number of grid points can be used in the integration. 
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FIG. 1. A schematic picture of the integration contour (dashed line) in the complex plane. The roots f2„ and poles of 
F(z) are marked with crosses. 



The RPA correction to the angular momentum operator p 



r(2) 



ij,kml 



)b\ l b J i, 



(14) 



(2) 

gives us the RPA contribution < J\ (f2„) > of the RPA solution f2„, where j\ are the single particle matrix elements 
of the angular momentum components i = x,y,z and the term b a bji is the boson image of the two quasiparticle 

(2) 

scattering contribution of the angular momentum operator. By choosing the weight function g(z) =< J> (z) > in 
the integral, ( pd| ) the corresponding RPA contribution to Ji can be calculated too. 

It should be pointed out that eq. (Q) takes into account also the contribution from possible spurious modes. In the 
paired case A > the function F(z) has a root at z = 0. This does not cause a problem for the integration since the 
factor z in the integrand in eq. ( |12| ) cancels this pole. After the transition to A = the spurious mode disappears. 

The above described method can be straightforwardly applied to obtain the correlation energy of the K — ground 
state band. The RPA for the excited configurations of the high- if bands needs special care concerning the selection of 
the "physical" phonon modes from the full RPA spectrum ±f2„ which implies positive and negative frequencies. The 
general criterion is to select from the pair of roots ±f2„ the one with the positive normalization. Such an identification 
is trivial for the ground configuration because there all Vl v > have positive normalization. However, for excited 
configurations the normalization check has to be done explicitly. Since the excited configurations under study are still 
close to the ground configuration the search for the possible physical phonons with f2„ < can be confined to a narrow 
energy region of a few MeV above zero. For this region the sum (|?]) is directly evaluated whereas for the higher lying 
part the loop integration is done. We mention that the correlation energy is also for excited configurations a negative 
(gain) term the value of which depends on the rotational frequency as well as on the configuration under study. 

In calculating the RPA contribution to the angular momentum components by means of eq. (|l4|), we observed that 
for some phonons, fi„, the backward going amplitudes become so large that the underlying quasiboson approximation 
(QBA) of the RPA is no longer valid. This happens in excited configurations near instabilities of the static pair 
field, which occur when two configurations with different A cross each other. Examples are the bands K = 7, 15 
(neutron pairing) and K = 15,22 (proton pairing), which will be discussed in section IVB. For the low lying RPA 
solutions we check explicitly whether the QBA condition, (y \ bybji \v) <C 1, is satisfied or not. We exclude roots with 
{y \ blfiji \ v) > 0.1 from the sum (|lj 



This procedure may be justified as follows. In the extended boson approximation (EBA) [[13 14 the non zero term 
{v\ b\jbji \ v) is approximatively taken into account by an iterative process. The relevant effect of the EBA as compared 



3 



to QBA is a strong reduction of the backward going amplitudes [jlj]. A full EBA treatment is too demanding for 
a realistic calculation like ours. Instead we use the rough approximation that the backward going amplitudes are 
quenched in the cases when QBA is not valid, i.e. if {v\ b'-ibji \v) > 0.1 the roots are removed from the sum @ of the 
RPA correction. These roots would give a very large (several units of K) contribution to the angular momentum. If 
they are included, Jrpa changes very rapidly near the crossings between configurations with different A. It is known 
the cranking model is unreliable near such crossings, where it seems justified to resort to an only rough correction of 
the QBA, which results in a smooth reasonable function Jrpa(<-<j)- 



C. Particle number projection 

The above treatment of the dynamical pair correlations by combining the HFB and RPA methods does not work 
for A — > 0. Another approximate method to include the dynamical pair correlations is the particle number projection 
(PNP) H||. In the PNP approach a variational state \N > with good particle number TV is formed by applying the 
projection integral 

\N> oc / d^e-^-^l > (15) 



to a TAC eigenstatc | > of the quasiparticle Routhian h' , eq. (^). Number projection accounts for gauge angle 
fluctuations of the pair field on a circle in the complex plane (Ae _l< ^) with the radius A, but disregards radial 
fluctuations. 

The PNP energy function to be minimized is the expectation value 

R PNP {A) =< N\H'\N > +\N (16) 
calculated with the projected state \N >, eq. ([l5|). We do not perform a full variation of the HFB amplitudes of the 
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reference state | > but minimize the energy (Hq) only with respect to the variable gap parameter A. The minimization 



of eq. (16) in PNP replaces the self-consistency condition (g) for calculating A in the HFB approach. Compared to 
the self-consistent pair gap (denoted by A HFB in the following) of an unprojected quasiparticle state, the pair gap 
A PArp optimizing the projected Routhian (|l^) is generally larger and stays non-zero through the phase transition 
region. 

In the HFB variant of TAC, the constraint < N >= N adds an additional dimension to the system of nonlinear 
equations given by eq. (|^) and the minimization with respect to e v and Performing the PNP the above condition 
is automatically satisfied. Instead A becomes another variable parameter of the Ritz variational problem, which is 
to be determined by minimizing the Routhian. However, the minimum practically coincides with the solution of 
eq. (|]), as illustrated in Fig. B This has the advantage that in the vicinity of the minimum the energy does not 
change very much. Thus the errors remain small if the minimum in A is not exactly found (or eq. (|]) is not exactly 
solved). Keeping A fixed greatly simplifies the actual calculation. However, one has to be careful that fixing A does 
not affect the general properties of the configuration too much. As will be discussed in Appendix [A|, neglecting the 
exchange term of the interaction may lead to problems in calculating the projected energy. Unlike the unprojected 
HFB function | >, the p rojected wave function (|l5|) is not be stationary. This may lead to problems that will be 



discussed in section II D 
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FIG. 2. Particle number expectation value in the unprojected case, < N >, and the total energy in the projected case, 
Epmp, as functions of A. This example corresponds to the yrast neutron configuration at u = 0.4 (MeV). The correct particle 
number < N > is represented by a horizontal dot-dashed line. 



D. Energy minimum in i? 

Since the shell correction part in total Routhian R, eq. (||), does not depend on the tilt angle z? one has 

dR d /dh'\ ,,, d , 

= i^^ / > = <^r>+2<|^— 1>. (17) 



dd d$ \dti/ 1 dti 

The term < > vanishes because | > is a stationary eigenstate to bl . Hence, the derivative of the Routhian 

becomes 

dR 1 3h! \ 

~M = \~M I = ~ w ( cos?9 < J i > - sirl tf < J 3 >) = -wJ±. (18) 
Thus. ! he requirement ^ = for the minimization is cast into the condition J± = 0. In other words, lu and J must 



be parallel 



ISp is not strictly valid, 
re condition for uniform 



In case of particle number projection, the state \N > is not an eigenstate to h! and eq. ( 
The energy minimum will no longer exactly agree with the condition of parallelity, which is t 
rotation. In the PNP calculations the minimum condition §f = needs to be fulfilled even if J± = is not satisfied. 

We generally found small differences (0° — 10°) between the values of § obtained from the condition J±(ll>) = 
and the energy minimum. Substantial deviations appear in regions of band crossings, where the cranking model is 
unreliable Jl5| anyway, and sometimes close to the band head. 



III. DETAILS OF THE CALCULATIONS 

We use the modified oscillator model with deformations S2, £4 and standard Nilsson parameters (l6j. The high- if 
bands in nucleus 178 W are found to be axially symmetric configurations. The tilt angle d is the angle between the 
symmetry (3-)axis and the rotational axis. The Strutinsky part E stru t(£v) of the TRS, eq. (^|) is obtained practically 
by including 8 oscillator shells for protons and neutrons, respectively. For calculating the quasiparticle term Rtac, 
eq. (H), it is sufficient to consider only a few shells around the Fermi surface. We included in the diagonalization of the 
quasiparticle Routhian the four N-shells closest to the Fermi surface. It is important to follow a certain quasiparticle 
configuration when seeking the minimum of R(u>, £2 £4, A, A, 1?). This is achieved by 'diabatic tracing': When changing 
one of the parameters determining h! the overlap of the quasiparticle wave functions with the ones before the step 
is calculated. By looking for the maximal overlap a one-to-one correspondence between the quasiparticle states is 
established. Figure || shows a quasiparticle diagram constructed in this way by using the step size Aw = 0.05 MeV, 
A^ = 5°. Keeping the occupation of such diabatic quasiparticle trajectories one usually follows a quasiparticle 
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configuration of a given structure. Problems may appear near quasi-crossings like the ones in Fig. 
configurations are mixed up. 



where the 
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FIG. 3. A quasiparticle diagram. The lines are connected by diabatic tracing, i.e. by finding the largest overlap between 
the points of two successive values of w. Solid lines correspond to quasiparticles with positive parity while the dot-dashed ones 
have a negative parity. Two quasi crossings are marked with double arrows. 

We have used two different strategies for finding the minimum with respect to A and $ for different w. 



1. 



Construct the configuration for A = 0. Choose a reasonable start value ds for different u>. 
Find minimum in § by tracing diabatically for fixed A. 
Find minimum in A by tracing diabatically for fixed d. 
Iterate by returning to (b) until total minimum is found. 

Construct the configuration for the band head for a guessed A. 
Find minimum for A. 

Go to the next w-point with fix A, tracing diabatically. 
Vary A and d as in lb - Id. 
Iterate by returning to (c) 



Usually 1. is better to work with because the minimum in d does not change much when pairing is added. Strategy 
2. is more appropriate for the ground configuration, because it has no counterpart in the case of zero pairing. 

The equilibrium shape of 178 W turned out to be rather stable. The deformation differs by less then 5% between 
different band heads (cf. |]l7| ). We also found that the deformation changes very little within the bands. Thus, we 
adopted the fixed deformations £2 = 0.229 and £4 = 0.034 calculated for the ground state. These values were used 
for all bands except K n = 25 + . This band contains a h 9 / 2 aligned proton, which drives the equilibrium values to 
£2 = 0.255 and £4 = 0.038. The relative energy of the bands is sensitive to the strength of the pairing force G 
In our calculations, G was fixed to match the values of the even-odd mass differences as calculated with PNP 



A„ = 1.15 MeV -► G v = 0.093 MeV 



A T = 1.23 MeV -> = 0.121 MeV 

In the PNP calculation the value of A is determined for each configuration at uj = 0.4 (MeV) by minimizing R(\) 
and kept constant for each configuration. In RPA and HFB calculations, A is adjusted according to eq. (Q) for each 
configuration and u>, because the HFB energy is much more sensitive to this parameter. 



IV. RESULTS 



The configurations of the considered if-bands are listed in table |jj They are in accordance with Q . The following 
presentation is divided into two parts. The first part treats the cases where HFB gap A HFB has collapsed already at 
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the band head (K^ = 25+ and 30 + ) while second comprises the cases where the pair field has transitional character 
(K* = 7-, 15+ and 22"). 

For a given configuration, there is usually only a weak systematic reduction of A when ui increases. It only changes 
substantially when a pair of protons or neutrons becomes aligned at a band crossing in a configuration with a low 
number of excited quasiparticles. For example, in the yrast band there is a drastic decrease of A„ through the crossing 
region at lu c « 0.3. It is well known! 15 that the cranking model has problems close to a band crossing. In this region, 
the different strategies (cf . section III ) of optimizing the mean-field parameters can result in different configurations 
and, as a consequence, in different self-consistent values of the parameters. In our TAC calculations for i^-bands, we 
observed the general tendency that the tilt angle $ increases monotonically from 0° at the band head to a value close 
to •& = 90° at highest frequency. In some cases there is a decrease in 1? at high w(> 0.5), which can be related to the 
crossing with other configurations. 





Neutron configuration 


Proton configuration 


T 


^{7/2+, 7/2"} 




15+ 


^{7/2+, 7/2-} 


^{7/2+, 9/2-} 


22~ 


^{5/2", 7/2+, 7/2", 9/2+} 


tt{7/2+,9/2-} 


25+ 


i/{5/2-,7/2+ 7/2-,9/2+} 


tt{1/2-,5/2+,7/2+,9/2-} 


30+ 


45/2", 7/2+, 7/2", 9/2+} 


tt{5/2+, 7/2+, 9/2", 11/2-} 



TABLE I. Occupied orbitals (fi lr ) of the multi quasiparticle excitations {IC). The quasiparticle configura- 
tions are taken from Q. Complete assignment: Neutrons 5/2" [512], 7/2+ [633], 7/2" [514], 9/2+ [624] and protons 
1/2" [541] , 5/2+ [402] , 7/2+ [404] , 9/2" [514] , 1 1/2" [505] 
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A. = 25+ and 30+ 



In these eight-quasiparticle configurations the static pair field has disappeared (A HFB = 0). The pair correlations 
are only of dynamical nature. One can expect the RPA to do a good job in describing them because one is sufficiently 
away from the critical region. As mentioned above, the PNP calculation leads to a finite value of A PNP which is 30% 
smaller then the values in the ground state (c.f. Fig. |I] and |^). There is a weak reduction towards higher u> within 
each band. 
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FIG. 4. A PArp as a function of u> for neutrons in the K n = 25 + and 30 + bands. Note, in these cases is A HFB =0. 
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FIG. 5. A PNP as a function of u for protons in the = 25 + and 30 + bands. Note, in these cases is A HFB —0. 



1. Routhians 

Experimental Routhians are calculated by means of the standard expressions given e.g. in |is|| . In order to remove 
the steep decrease of R(oj) a term ZQMeV^uo 2 is added to all R values, which makes the differences between the 
various curves better visible. The calculated ground state energy R g . s .(u) = 0) is set to zero for all calculations and 
the experiment. The PNP and RPA calculations give an energy gain relative to the HFB results. Due to the blocking 
effect, this gain is normally smaller in the excited states than in the ground state. As a consequence the energy 
relative to the ground state gets larger for RPA and PNP than for HFB calculations. This systematic tendency is 
substantiated by the examples to be shown. 

For the = 25 + band the RPA yields a substantially better agreement with experiment than PNP (cf. Fig. |^). 
The larger discrepancy for the Jf" = 30 + may reflect some inaccuracy of single particle levels. Although the HFB 
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calculation gives a good estimate of the angular momentum [g] it underestimate the relative energy of the bands. The 
rotational frequency of the band head is also quite well reproduced by RPA and PNP. 
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FIG. 6. The energies in the rotating frame R(u>) relative to the ground state for the JC = 25 + (upper) and JC = 30 + 
(lower) bands. A term of 30a; 2 is added to the R values. The experimental values are taken from Q. 

In Fig. [t] the energy differences between the paired (PNP or RPA) and unpaired calculation (Hartree-Fock) are 
shown. They characterize the effect of the pair correlation onto the energy. We can see that the RPA gives a larger 
pair correlation energy then PNP in accordance with ^0|. This indicates that the RPA is a better approximation 
than PNP in the A HFB = regime. One notices that the correlation energy is only weakly reduced by rotation. 
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FIG. 7. The total pair correlation energies in the rotating frame Rp a ired(w) ~ Runpaired for the the K n = 25 + and 30 4 







2. Angular momentum 



The investigation demonstrated that the angular momentum and the dynamical moment of inertia of 178 W can 
be understood by assuming that the nucleons move in a rotating mean field with no pairing. The strong reduction 
of the mome nt of inertia as compared to the rigid bodyf| value is due to the nuclear shell structure (cf. Q and 



IV C). As seen in Fig. H, the inclusion of dynamic pair correlations does not change this result. Note that a 
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linear term has been subtracted so that the differences between the curves are considerably magnified. 
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FIG. 8. Angular momentum J(uS) for K w = 25 + and 30 + . A term of 50a; is subtracted from the J values. The experimental 
values are taken from Q. The RPA and HFB results practically coincide in this picture. 

All approaches reproduce the experimental dynamical moment of inertia {J^ — g~) which is substantially below 
the rigid body value. The differences of the function J(to) between paired and the unpaired calculations are only 
marginal. The dynamical pair correlations do not change the angular momentum very much, but they do increase 
the energy difference between bands with a different number of blocked quasiparticles. 



3. Branching ratios 



In the TAC model, the B(M1) and B(E2) values are calculated by means of the expressions 
B(M1) = — [sintf (J 37T + 2.35SW - 2.24S 3 u) - cost? (J l7r + 2.35SW - 2.2AS lu )] 2 

87T 



(19) 



1 This limit is referred to the rotational spectrum E(I) = A r i g I(I + 1) — K 2 
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and 
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(sintf) 4 ^ 



(20) 



where J, S and Qo are the expectation values of the angular momentum, the spin and the quadrupole moment, 
respectively, as calculated with the TAC states. It is common practice to present the experimental branching ratios in 



the form 



(gK-gn) 



which is obtained assuming that the strong coupling limit plfl is valid. We choose to display 

the calculated branching ratios (which of course do not rely on the strong coupling assumption) in the same way. The 
theoretical ratios are obtained as 



(gi< - 9r) 
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(21) 
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where K, the value of the angular momentum at the band head, is kept constant and J is the calculated value of the 
angular momentum. Qo is chosen as in 
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FIG. 9. The ratio 
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I q 9 for the K 7 ' — 25 + and 30 + bands. The experimental values are taken from 



Observe the 



The theoretical and experimental ratios are given in Fig. |9| Since the ratio depends strongly on the orientation of 
the proton and neutron angular momentum the good agreement confirms our calculated geometry. The experimental 
errors of the branching ratios are too large to discriminate between the PNP and HFB calculation. 



B. K n = 7", 15+ and 22" 

The if" = 7~, 15 + and 22~ bands are analyzed and presented in the same way as the = 25 + and 30 + bands. 
These bands, which correspond to the excitation of two to six quasiparticles from the K = ground configuration, 
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have a reduced but non-zero static HFB pair field at the band head. The PNP is expected to be stable while the RPA 
can run into problems since one is close to the transition A HFB — > 0. 
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FIG. 10. A HFS (grey) and A PNP (black) as a function of a> for neutrons in the K n = 7", 15+ and 22" bands. 
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FIG. 11. A HFS (grey) and A PNP (black) as a function of to for protons in the K 71 = 7 _ , 15+ and 22" bands. 

Figure |o| and [ll] show the calculated values of A. The values of A PNP are similar to the ones for the K - * = 25 + 
and 30 + bands in Figs. ^ and ^, except K w — 7", for which the protons remain in the ground configuration. For the 
two quasi neutron bands K n = 7" and 15+ the gap A^ FB is reduced strongly already at the band head and becomes 
zero at u> ~ 0.3 (M eV). T hese values indicate that the pair filed is nearly instable, which causes problems to be be 
discussed in section 1VB 2. Though A BFB is strongly reduced in the two quasi proton configurations = 15 + and 



22 it changes only weakly with u>. 



1. Routhians 



Both the PNP and RPA calculations fairly well reproduce the Routhians at high u>. There, again the tendency 
is seen that the dynamical pair correlations enlarge the energy distance between the bands. The exception is the 
= 15 + for uj > 0.38 MeV. This discrepancy will be discussed in section IV B 2. 
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FIG. 12. The energies in the rotating frame R(uj) relative to the ground state for the K n — 7 (upper), 15 + (middle) and 
22~ (lower) bands. A term of 30cj 2 is added to the R values. The experimental values are taken from [?]]. 

Figure [l3] shows that the pair correlation energies for these lower K values are larger and depend stronger on the 
rotational velocity than the ones in Fig. fjj The HFB + RPA calculations cannot reproduce the backbend in the 
K"* = 7~ and 15 + bands since A BFB = for uj > 0.25 MeV. It is seen as the kink in Fig. [f| and the upbend of 
the PNP calculation in Fig. [l^. The erratic behavior of the PNP pairing energy at large ui can be understood by the 
fact that the energy surface has a shallow minimum and is therefore easily disturbed by structure effects caused by 
quasiparticle crossings. Most of these crossings do not exist in the HFB + RPA calculations because A^ FB = 0. As 
a consequence, R depends more smoothly on uj, as the experiment does. We consider the fluctuations of R as an error 
of the PNP approach, which mimics the dynamical correlations by a static A PNP . 
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FIG. 13. The total pair correlation energies in the rotating frame R P aired(u) — Runpavred f° r the the K" = 7 , 15+ and 22 . 

2. Angular momentum 

As seen in Fig. |l4[ the calculations give reasonable agreement with the experiment. The neutron backbends in the 
K 71 = 7~ and 15 + are not reproduced in the HFB + RPA calculation because the static pair gap collapses already 
before the backbend. This feature could be corrected by using a somewhat larger pairing strength G v . The PNP 
calculation gives approximately the right value of lo c for the neutron backbend but gives a slightly too low value of 
J(u) for the K n = 7~ band. 

In the K 77 = 15 + band, all the calculations underestimate the angular momentum by ~ hTi after the backbend at 
u> ~ 0.38 MeV. We assign this increase to a configuration change caused by a crossing between a hg/2 proton orbital 
and a hn/2 proton orbital. As a result a K 71 — 19 + band continues the K* = 15 + band. This crossing occurs at 
the same u> as the neutron crossing causing the normal backbend. The same proton crossing occurs in the K v = 22~ 
band at u> ~ 0.48 MeV which continues as a = 26~ band. It is seen in Fig. |lj as an upbend. The configurations 
of the K 71 = 19 + and 26 _ bands are listed in table ||. The calculations for the if" = 7~ band after the backbend also 
gives a bit too low angular momentum. It cannot be explained by the proton crossing discussed. This may indicate 
that the missing angular momentum in the K n = 15 + and 22~ bands has another origin than suggested above. 
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Neutron configuration 


Proton configuration 




19+ 


^{7/2+ , 7/2"} 


^{1/2", 7/2", 7/2+, 9/2" 
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26- 


t/{5/2-,7/2+,7/2-,9/2+} 


7r{l/2-,7/2-,7/2+,9/2- 


} 
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TABLE II. The quasiparticle configurations used in the article are taken from J7|. (Neutrons: 
5/2" [512], 7/2+ [633], 7/2" [514], 9/2+ [624] Protons: 1/2" [541], 7/2" [523], 7/2+ [404], 9/2" [514]) 
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3. Branching ratios 



The good agreement with experiment in Fig. 15 for the K 7 ' = 15 + and 22 confirms the geometry of our calculations. 



Again, the large experimental errors of the branching ratios do not allow to distinguish between the PNP and HFB 
calculations. 
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In the case of the K"" = 7~ band, the calculation with A = reproduces the branching ratios better. We consider 
this as accidental, because the presence of the upbend seen in the I{uS) curve of Fig. |l4| is a clear indication of a 
finite pairing gap for the neutrons. The protons are definitely in their ground configuration with a non-zero gap. 
The K v — 7~ band has small B(M1) values compared to other JC-bands. This is expected because the normally 
dominating proton contributions to the B(M\) value ( |T9| ) are small for the ground configuration and very sensitive 
to the size of A^-, It is largely canceled by the negative neutron spin contribution. Hence, the resulting ratio is not 
reliable. 



C. COMPARISON OF ALL If -BANDS 

In Figs. |l6|-|l8| we compare the Routhians and the angular momenta of different if-bands at a fixed frequency 
u) = 0.45 MeV. The K — 30 band is left out from the plots since it starts at a higher uj. The bands with different 
orientation of J have approximately the same energy R in the rotating frame. Their energies in the laboratory frame 
for a given value of J do not differ much as well. That is so, because the energy needed for generating angular 
momentum along the two principal axes 1 and 3 is nearly the same. This feature, which is reflected by many bands 
with different orientation being close to the yrast line, is characteristic for the (N,Z) region where 178 W is situated. 
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FIG. 18. The angular momentum component along the 1- and 3-axis for different K at u = 0.45 as calculated with HFB. 
The K- values plotted are O(s-band), 7, 15, 22 and 25. 
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FIG. 19. The negative parity states in the 50-82 protons shell. The larger gaps between the levels in the upper part of the 
shell promote deformation alignment of angular momentum. The levels are labeled with the asymptotic quantum numbers 
[Nn z mi]Q. 

In the upper and middle part of the 50-82 proton and 82-124 neutron shells there are many orbitals with higher 
fl values that strongly couple to the prolate deformed potential whereas the low-f2 orbitals are occupied. Figure |l9| 
shows the /in/2 orbitals as an example. For the high-f2 orbitals the Coriolis coupling between states with and fi + 1, 
which generate a contribution to Ji, is weak, because the energy distance is large and the coupling matrix-element 
(~ y/ j 2 — fl 2 ) is small. This is very different from the nuclei situated at the bottom of the shell. For the low- orbitals 
the energy distance is small and the coupling matrix element large. Hence it is easy to generate a contribution to J%. 
Thus, the strong coupling suppresses the angular momentum alignment along the collective 1-axis, which causes the 
reduction of the moment of inertia relative to the rigid body value. 

On the other hand, it is energetically favorable to generate a contribution to J3, because there are many orbitals 
with a large |f2| near the Fermi surface and it does not cost much energy to generate a particle-hole excitation with 
a large K = il p — (— fi/j). At the bottom of the shell these excitations cost much more energy, because the orbitals 
with large |0| are far away. The net result is that at the beginning of the rare earth region the moment of inertia is 
larger than the rigid body value and generating angular momenta along the 3 axis is unfavorable. In the upper part 
of the region, the moment of inertia is smaller than the rigid body value and generating angular momenta along the 
1 and 3 axis is equally favorable (see and |2^]). The inclusion of pairing does not change this conclusion. 

Figure |l8| shows the J\ and J3 components for different i-T-bands calculated with HFB at a typical frequency 
uj = 0.45. Within ±2h, the component J3 is equal to K which means that all bands are not far from the strong 
coupling regime. The largest deviation is found in the K — 15 band where J3 ss 17 at uj = 0.45. This is an effect of 
the configuration mixing in the neutron system, which is induced by quasi crossings in a region of high quasiparticle 
level density. The total spin J is almost independent of K, in accordance with the discussion above. Generating the 
components J\ and J3 of angular momentum costs a comparable amount of energy. The slightly different behavior of 
the K = 25 band is due to extra rotational alignment from the aligning hg / 2 proton orbital, see discussion in Q . 

Summarizing our calculations concerning an improved treatment of pairing we find that i) the conclusions of earlier 
work Q can be confirmed and ii) the angular momentum J{uj) is insensitive to the pair correlations in the frequency 
region after the first band crossing. This is at variance with the results of ||[| for low-K bands, where a reduction 
of several units of angular momentum by the dynamical pair correlations was found. The reason is that the single 
particle orbits near the Fermi surface, which contribute most to the collective angular momentum, are blocked for the 
pair correlation by generating the high K. The ratio \^^^-{oj)\ is also reasonably well reproduced which support 
the calculated tilted geometry of our solutions. 



We included dynamical pair correlations into the tilted axis cranking approach in two different ways: The pairing- 
RPA method, which allows for harmonic vibrations on top of the HFB mean field, and the particle number projection 
which describes the dynamical correlations by an increased static pair gap in conjunction with projection onto good 
particle number. We studied the high iC-bands in 178 IF by means of both methods. As known from investigation of 
the band head energies E), the HFB approach tends to underestimate the energetic separation between the bands. 



V. CONCLUSIONS 
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Inclusion of the dynamic correlations improves the energy of bands with zero static pair correlations relative to each 
other and to the bands with pronounced static pair correlations. When the static pair field is substantial we could 
not find a systematic improvement of the relative energy due to the dynamical pair correlation. This seems to be in 
contrast to the result for K = bands in H||. 

The pronounced reduction of the moment of inertia of high- if bands relative to the rigid body value, seen in the 
experiment, is not due to the pair correlations. It is caused by orbitals close to the Fermi surface which in the mid 
and upper part of its shell is disfavoring the generation of angular momenta perpendicular to the symmetry axis. It is 
typical for this region, that the generation of angular momentum along the symmetry axis by particle-hole excitation 
is equally favorable as compared to the rotational alignment i.e. collective rotation along the perpendicular axis 
(see |p[ and ^2j). The inclusion pair correlations, both static and dynamic, does not significantly affect this feature 
of high-if bands, which is in stark contrast to low- if bands. 

In bands where the static pair gap of the HFB treatment has collapsed to zero the RPA method is simpler to use 
than PNP because there is no need to minimize the Routhian with respect to the A and A parameters. The pair gap 
introduced in these PNP calculations causes irregularities in the energy and angular momentum which are caused 
by quasiparticle band crossings. The RPA method, which does not have these crossings, describes the experimental 
results better. Thus, we consider the irregularities appearing in PNP as an artifact of the approximation and RPA 
superior in this region. On the other hand, when the HFB mean field pair gap is reduced but not yet zero the RPA 
method has problems because of its deficiencies near the transition to zero pairing. In these cases it turns out to 
be important to check the quality of the quasiboson approximation when calculating the RPA contribution to the 
angular momentum. The PNP method is more stable in these cases. 

We gratefully acknowledge valuable discussions with R. G. Nazmitdinov. Work is supported by the DOE grant 
DE-FG02-95ER40934. 
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APPENDIX A: THE ROLE OF THE EXCHANGE TERM IN PARTICLE NUMBER PROJECTION 



The HFB pairing energy is usually calculated in Hartree approximation, i.e. neglecting the exchange term by 
factorizing the pairing matrix element < P^ P > « < P^ >< P > . This approximation is justified for the calculation 
of the pairing energy contribution without PNP. However, it was recently suggested [|3| that the neglect of exchange 
terms in performing PNP can lead to dangerous poles in the resulting total Routhian surface (TRS) Such an unphysical 
behavior of the PES was indeed found in our calculations and it was traced back to the above mentioned factorization. 
For the sake of completeness we sketch the argumentation presented in more detail in E3] . 

and the canonical (BCS-like) form || 



Using in the construction of the PNP state 
>, the number projected pairing energy E pa i r 



N > eq. 

= -G < 
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of the HFB state 
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and the usually neglected exchange term 
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These expressions contain the canonical BCS-amplitudes Uk,Vk of the quasiparticlc state and the matrix elements 
Pkk' of the pair operator in the canonical basis. The bracket 



<ip = 0\ip>= n fe>0 (ill T < ; 



22i V 



i) i..e 



(A4) 



is the overlap function between gauge rotated quasiparticle states. 

One may encounter a zero denominator (u\ + v^e 2 ^) in both energy terms (A2 A3) e.g. for u% = vu and ip = -|. 
This is because the dou ble zero in the denominator cannot be canceled by the corresp onding single zero in the overlap 
< ip — 0\ip >, eq. (A4). However, when summing up the two contributions (A2 A3) to the full pairing energy E pa i r 
such unphysical poles do exactly cancel. 
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FIG. 20. The energies E pair and the Ep a % ect calculated for u = 0.3 in the K n = 7" band with X N = 49.54. Both curves are 
normalized to zero at A = 0.4. 



Our calculations confirm the conclusion that a reliable calculation of the TRS with PNP has to be done with the 
full expression. In Fig. |2^ the full neutron energy is sh own as a function of the neutron gap (A) and it is compared 
to the one where only the direct pairing energy term ( |A2| ) is taken into account. The full energy has the expected 
parabolic shape with the minimum whereas the curve of the direct term alone displays an unphysical pole around 
A w 0.6 MeV. 
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FIG. 21. A determined from E^tr* and E vair , respectively, for the K n = 7~ band, see also Fig 
small difference except at oj = 0.5 where there is a pole in the energy function Ep££, ct , see Fig. t 
found keeping all other parameters constant at the self-consistent values found using the full projection. 



| and [l^. There is only a 
These values of A were 
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FIG. 22. Routhians R calculated with £7*^°' and E vair , respectively, for the K* — 7~ band, see also Fig. [l2| There is only 
a marginal difference between both Routhians where all other parameters are kept constant at the self-consistent values found 
using the full projection. 

Such a strange behavior does not happen often and usually there is only a minor difference between the extracted 
A-value at the minimum for the full and direct energy (cf. Fig. |2l]). Observables like the energy are not strongly 
affected by the exchange term except close to poles (cf. Fig. ^2|) . With the full pairing energy, we obtain a minimum 
for a slightly different A and, therefore, a different G is needed to match the experimental A-values. We found that 
a Gjr of 0.121 MeV instead of 0.119 MeV should be used for the protons while the G„ did not change when using the 
full expression for the pairing energy. 

The probability of accidentally hitting a pole is not large but it happened a couple of times in our calculations. 
The tail (cf. Fig. |2^) of a pole can also affect the results and this is, of course, much harder to detect. The energy 
surface calculated with only the direct term jumps when passing through the pole. This is because a pole has gone in 
(or out) to the area in the complex plane around which we are integrating. The pole would turn into a step function 
if it was possible to perform the integration exactly. In order to avoid such unphysical one should generally apply the 
full expression of the PNP pairing energy. 
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